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The reaction 7 + p — > K + + E + tv was used to determine the invariant mass distributions or 
"line shapes" of the E + 7r~, E~7r + and E°7r° final states, from threshold at 1328 MeV/c 2 through 
the mass range of the A(1405) and the A(1520). The measurements were made with the CLAS 
system at Jefferson Lab using tagged real photons, for center-of-mass energies 1.95 < W < 2.85 
GeV. The three mass distributions differ strongly in the vicinity of the 1 = A(1405), indicating the 
presence of substantial 7 = 1 strength in the reaction. Background contributions to the data from 
the E°(1385) and from A"*E production were studied and shown to have negligible influence. To 
separate the isospin amplitudes, Breit-Wigner model fits were made that included channel-coupling 
distortions due to the NK threshold. A best fit to all the data was obtained after including a 
phenomenological 1 = 1, J P = l/2~ amplitude with a centroid at 1394 ± 20 MeV/c 2 and a second 
1=1 amplitude at 1413 ± 10 MeV/c 2 . The centroid of the 1 = A(1405) strength was found at 
the E-7T threshold, with the observed shape determined largely by channel-coupling, leading to an 
apparent overall peak near 1405 MeV/c 2 . 



PACS numbers: 13. 30. Eg 13.60.Rj 14.20.Gk 25.20.Lj 

I. INTRODUCTION 

The A(1405), situated just below NK threshold, has 
been an enigmatic state in the spectrum of strange 
baryons for decades. First seen in bubble chamber exper- 
iments in the 1960's [1], there have been remarkably few 
measurements of this state to date. The most prominent 
feature of the state is that its invariant mass spectrum, 
which we call the "line shape", has always been seen to 
be distorted from a Breit-Wigner form, indicating that 
there are strong dynamics at work that are not seen in 
more typical resonances. Almost all theories agree that 
this is due to the state's strong coupling to NK, but the 
exact nature of this coupling is as yet unknown. Due to 
its mass being below the NK threshold, it is not possi- 
ble to produce it directly in kaon beam experiments, so 
accessing this state experimentally has been a challenge 
compared to other strange baryon resonances. Precise 
measurements of the line shape should yield information 
on what dynamics play a significant role in the A(1405), 
and lead to a deeper understanding of the additional am- 
plitudes that may exist in this mass region. 

A. Theories of the A(1405) 

Explaining the mass of the A(1405) has also proved 
to be a problem. The state does not fit well within the 
constituent quark model that has otherwise worked re- 
markably well for understanding the masses of low-lying 
baryon resonances [2] . Theoretical investigations into the 
nature of the A(1405) were discussed from the days of its 
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prediction by Dalitz and others [3], and there has been 
a surge of interest in recent years. Chiral unitary the- 
ory [4-6] combines chiral dynamics with unitarity con- 
straints based on effective meson-baryon interactions. In 
this class of models, the A(1405) is dynamically generated 
as a rescattering of all pseudo-scalar meson and octet 
baryon states that couple to it. Definite predictions have 
been made of what the line shape of the A(1405) should 
be for photoproduction near threshold [7] . In this model 
the interference between a dominant isospin 1 = am- 
plitude and a smaller 1 = 1 amplitude causes the line 
shapes for each Y,ir channel to be different. 

Further developments of the chiral unitary approach 
have shown that the A(1405) may be composed of two 
1 = poles, whose couplings to various particle final 
states and whose initial-state populations differ according 
to the reaction under investigation [8-10]. The A(1405) 
plays a special role in these theories as the archetype 
of a dynamically generated rescattering state, but the 
models also impinge on the nature of non-strange nucleon 
resonances such as the iV(1535)Sii and A r (1440)Pn. 

In another approach [11] the pair of states £(1385) and 
A(1405) are treated together in a kaon double-pole model 
with an explicit assumption of photon dissociation into a 
real K + and a virtual K~ as illustrated in Fig. 1. This 
model also made specific predictions for the mass distri- 
butions of £-7r and A7r° final states that we will mention 
later. In a meson exchange model [12] the A(1405) is 
generated dynamically by the coherent addition of lj, p 
and scalar-meson exchanges. Here, too, it appears as a 
two-pole structure in the NK Sqi partial wave. 

Other theories see the A(1405) as a bound state of 
NK [13] alone, and the dynamics that create it are 
expected to have significant repercussions on whether 
bound K~pp states exist. In other views, the A(1405) 
is pictured as a true 3-quark state [14], or as a negative- 
parity [15] or positive-parity [16] hybrid state. 

The spectrum of / = 1 £ excited states is also pre- 
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dieted to be quite different in different models. As dis- 
cussed recently in Ref. [17], standard so-called quenched 

quark models put the J p = \ E at about 1650 MeV/c 2 , 
while an unquenched model that allows [qq] [qq] q S- 
wave configurations expects it near 1380 MeV/c 2 . Some 
evidence for a light negative-parity E has been dis- 
cussed for several years [18-21]. In a different vein, 
meson-baryon dynamical models place a E* (l/2)~ 
near 1430 MeV/c 2 [22], 1475 MeV/c 2 [23], or 1620 
MeV/c 2 [24]. It is therefore of interest to look for ev- 
idence of isospin one strength in the E-7T system in the 
same neighborhood as the isospin zero A(1405). Discov- 
ery of a resonant structure in that same mass range could 
be decisive in picking among models of baryonic excita- 
tions. 



B. Experiments on the A(1405) 

While there has been continual theoretical interest in 
the A(1405), there have been remarkably few measure- 
ments made of this state, known for more than half 
a century and given a 4-star rating by the PDG [25]. 
Bubble chamber experiments using hadronic beams at 
Brookhaven [26] and CERN [27] have long been the only 
experiments to identify the line shape, and that with 
barely adequate statistics. The mass and width of the 
A(1405) cited by the PDG are based primarily on these 
measurements. 

In recent years, with the development of higher- 
statistics experimental capabilities, there has been a re- 
newed interest in measuring the A(1405). These include 
measurement of the E°7r n line shape in proton-proton 
collisions at COSY [28], measurement of the line 
shapes in proton-proton collisions at HADES [29], and 
again the measurement of the E*^ line shapes using 
photoproduction by LEPS [30, 31]. The LEPS mea- 
surement has some overlap with the energy used in our 
measurement, but due to limited statistics, their re- 
sults depended on broad averaging over kinematics and 
line shape comparisons with existing theoretical curves, 
rather than on new and decisive fits to the data. 

In this paper, we report results of a measurement with 
large statistics accumulated with the CLAS system in 
Hall B of Jefferson Lab. With good mass resolution for 
the A(1405), we will show for the first time a measure- 
ment of all three E-7T line shapes. The center-of-mass en- 
ergies (W) in this experiment covered a wide range from 
near production threshold of the A(1405) up to 2.85 GeV, 
which allowed us to measure the energy-dependence of 
the line shapes. The results are shown after summing 
the line shapes over all kaon production angles for each 
energy; another paper that is in preparation will show 
the differential cross sections for each energy [32] . 

The line shapes are differential in the E-7T invariant 
mass, m, and extracted for nine bins in the initial state 
7P invariant energy W. We anticipate that the W depen- 
dence of the 1 = and 1 = 1 contributions varies slowly. 



This is because the data stem from an associated pro- 
duction experiment accompanied by a kaon, as opposed 
to a direct formation experiment, so that the connec- 
tion of energy W to the properties of the excited strange 
rcsonance(s) is indirect. This is indicated in Fig. 1, for 
example, where the intermediate hyperon is construed to 
be created via an off-shell kaon as part of a t-channel 
interaction between the incoming photon and the target 
proton. 




FIG. 1. (Color online) Creation of the three-body K + E 7r 
final state via an intermediate hyperon in the reaction 7 + 
p — > K + + E + 7r. In this particular example, a t-channel 
exchange enables an off-shell kaon to create a A(1405) that is 
sub-threshold for on-shell NK reactions. 

The outline of this paper is as follows. Sections II- 
IV present the details of the setup of the experiment, 
the event selection used to extract the yield of the vari- 
ous final states, the acceptance corrections and the data 
normalization. Section V describes our method of ex- 
tracting the yield of signal events from the data, Sec- 
tion VI presents the line shape results for the three E7r 
final states, and Section VII discusses the systematic un- 
certainties and the methods used to test the reliability 
of the measurements. Section VIII explains the method 
used to fit the line shapes we obtained. In Section IX the 
outcome of fitting the mass distributions with this model 
is given, and we conclude with Section X. 

II. EXPERIMENTAL SETUP 

The data for this experiment were obtained during 
May and June of 2004 with the CLAS detector, located 
in Hall B at the Thomas Jefferson National Accelerator 
Facility The run, known as glla, used a 40 cm unpo- 
larized liquid hydrogen (LH2) target and an incoming 
unpolarized real-photon beam. Bremsstrahlung photons 
with an endpoint energy of 4.019 GeV were created via 
the CEBAF accelerator electron beam and a 10 -4 radi- 
ation length gold foil. Electrons that radiated a photon 
were identified with the CLAS tagger [33] to obtain en- 
ergy and timing information between 20% and 95% of 
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the endpoint energy. 

Details of the CLAS detector can be found in Ref. [34]. 
Here we give a very brief description of the main compo- 
nents used in our analysis. CLAS was equipped with a 
superconducting toroidal magnet with six identical sec- 
tors surrounding the beam line. The field was selected 
to bend positive particles away from the beam line. A 
34-layer drift chamber system in each sector provided 
charge and momentum information for charged parti- 
cles. Momentum resolution Sp/p was « 0.5%. The target 
was surrounded by a 24 element plastic scintillator Start 
Counter used in the trigger to select charged tracks leav- 
ing the target. Finally, a system of 342 time of flight 
(TOF) scintillators was used in the trigger which also 
determined the duration of flight of each charged parti- 
cle. For the glla run period, the trigger required a hit in 
the tagger system in coincidence with Start Counter and 
TOF hits in at least two of the six sectors. A sector trig- 
ger required hits in a Start Counter paddle and a TOF 
paddle within 150 ns of each other. With this setup, the 
glla run accumulated over 20 billion events, including a 
large sample of excited hyperon states. More details of 
the setup and analysis can be found in Ref. [35] . 



III. EVENT SELECTION 

For reference, we include Fig. 2 to illustrate the various 
data-handling paths in this analysis. The main interest 
lies in the T,ir hnal states of the hyperon decays, par- 
ticularly A(1405) — > YiTi. The analysis channels can be 
divided into two main categories depending on the final 
charged particles detected. The first case is when a K + , 
p, and 7T _ are detected, while the second is when a K + , 
7r + , and 7r~ are detected. The main strong final states 
of interest are then A+E+Tr", A+£°7r°, K + T,-it+, and 
A'+Att . The latter is mainly due to £(1385) -> Att°, 
which is one of the significant backgrounds for isolat- 
ing the A(1405). Since the A and £ hyperons decay 
via the weak force, we will call the final states including 
a hyperon the strong final states. During analysis the 
hyperons were reconstructed through their weak decay 
products. The ground-state hyperon decays £ + — > pir°, 
£ + — » nir + , £° — > 7 A — > jpTT~ , £~ — > nir~ , and 
A — > pn~ are detected in our analysis. 

The analysis procedure was as follows. The raw set 
of reconstructed and calibrated events was the same as 
used in several previous CLAS publications [36-39]. We 
selected events with all of the required charged particles 
for each channel of interest. After some cuts to reduce 
backgrounds, described below, a kinematic fit was ap- 
plied when all final state particles but one were detected. 
Otherwise, the missing mass squared was computed for 
the case of the £°7r° channel. From this, the ground state 
hyperons of A or £ were reconstructed and selected for 
analysis in each channel. For brevity, we will label the 
channels of interest by their ground state hyperon-pion 
combination, i.e., £+7r _ , £+7r _ , £°7r°, £ - 7r + , and A-7T , 



where the two £ + decay channels are distinguished with 
a subscript denoting the final-state baryon. Analyzing 
both £ + channels served as a cross check of uncertain- 
tics introduced by our analysis methods. 

A. Initial selection of particles 

The effects of selection cuts discussed below are sum- 
marized in Table I. In this analysis all channels of interest 
have a final state K + . For the entire data set, the masses 
of charged particles were calculated from the momenta 
given by the drift chamber tracks, and the timing given 
by the RF-corrected tagger timing and TOF scintillator 
timing. A loose mass cut was made to select events with 
a kaon candidate, and these were retained for analysis. 
Events were then required to have all charged particles 
reconstructed in the fiducial region of the detector, and a 
few malfunctioning timing detectors were identified and 
removed. The fiducial region of good acceptance and 
good Monte Carlo matching was the same as used in 
previously-published analyses [36-39] of the same data 
set. Copious non-strangeness events were removed by 
testing the hypothesis that a candidate K + was actu- 
ally a 7r + or a proton. The leading backgrounds were 
7P — > pK + ir~ with nothing missing, jp — > pir + ir~ (n ) 
with the 7T° missing, and jp — > 7r + 7r + 7r~(n) with the 
neutron missing. 

Previously-established corrections to the reconstructed 
tracks were applied, such as momentum corrections for 
small imperfections in the magnetic field map, and energy 
losses due to the charged particles traveling through the 
target and detector material. The incident photon en- 
ergy was corrected for the known mechanical sagging of 
the tagger hodoscope. In all cases a timing cut was ap- 
plied to remove events where a tt + was misidentified as 
the A + (see below). A primary event vertex cut along the 
beam direction selected events very cleanly from the LH2 
target and rejected events from foils. Minimum momen- 
tum cuts based on the identity of the particle (0.3 GeV/c 
for protons and K + , 0.1 GeV/c for n ) were applied. 

As most of the background in these channels came from 
strangeness-free events, cuts on the timing of particles 
were crucial to correctly select kaons. In the CLAS de- 
tector, the distance a charged particle travels through 
the drift chambers (I), the accelerator RF- and vertex- 
corrected event start time (to), and the time that the 
particle hit the TOF paddles (ii) were recorded along 
with the particle's magnitude of momentum (p). From 
this information, the measured travel time was calculated 
as 

^nieas — ^1 — ^0- (1) 

Alternatively, we assumed a mass hypothesis for the par- 
ticle, mo, and used the measured momentum to calculate 
the velocity of the particle as 

/?calc = / „ ^; (2) 

VP + m o 
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FIG. 2. (Color online) Representation of reaction and decay channels used in this analysis. Particles X given in the columns 
were reconstructed via kinematic fitting or the missing mass determination. Branching fractions are indicated as percentages. 
The right-most column lists the figures related to the given channel. 



TABLE I. The number of events remaining after each selection cut, in 1000's. 



Selection 


Channel 

E + 7T-, ATT , S°7r° E + 7T", T,-tt + 


Detected particles 


K + pir 




Initial kaon selection 


64,026 


35,627 


Fiducial cuts 


31,486 


16,662 


Remove false K + due to ix + or p 


4,852 


10,045 


Loose ATOF cuts 


3,093 


6,576 


Vertex z cut 


3,066 


6,464 


Minimum \p\ requirements 


3,047 


6,233 


Precise ATOF cuts 


2,415 


3,912 


Kinematic fit or MM 2 cut 


818 


233 


1,052 


Selection on ground state hyperon 


440 238 


76 


316 338 



and together with the reconstructed flight distance I de- 
termined the calculated flight time of the particle as 



B. Selection of events for analysis 



^calc 



/ 



(3) 



Taking the difference between these two timing measures 
gives 



ATOF = t n 



^calc ■> 



(4) 



and cuts were applied on this quantity as a function 
of particle momentum. Figures 3 and 4 show the 
momentum-dependent cuts applied to select the tt~ , 7r + , 
and K+. 



In all channels, the data were divided into 10 bins 
of energy spanning 100 MeV in the center-of-mass en- 
ergy W (= y/s), and 20 angle bins in the center-of- 
mass kaon angle. All selection cuts and fits to the data 
were done independently for each bin. For channels 
with the J^ ± TT zfl combinations, kinematic fits were applied 
with fixed mass of the undetected neutron or ir° (one- 
constraint or 1C fits). This optimized the information 
based on the measured momenta while balancing the en- 
ergy and momentum of the reaction. 
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FIG. 3. (Color online) Difference in particle time-of- flight, 
ATOF, versus the measured magnitude of momentum for the 
7r + (top) and 7r~ (bottom) for a given energy bin of 2.35 < 
W < 2.45 GeV. The different horizontal bands correspond to 
the 2 ns time structure of the CEBAF beam. The magenta 
lines show where the cuts were applied to select each particle. 




Ipl (GeV/c) 



FIG. 4. (Color online) ATOF versus the measured magnitude 
of momentum for the K + candidates at several analysis stages 
for a given energy bin of 2.35 < W < 2.45 GeV. From top to 
bottom, the plots correspond to: after the K + misidentifica- 
tion rejection cut; after selecting the detected in-time n~ as 
shown in Fig. 3; and selection on both the detected n~ and 
7r + . The magenta lines represent the selection cut on the K + . 
The last spectrum is seen to be much cleaner after the ATOF 
selection cuts on the pions. 



1. Event selection for Kit and Ep7r 

In these channels we reconstructed the particles K + , 
p, and 7T _ , with a missing ir°. The 1-C kinematic fit was 
applied to the selected particles, and those events with a 
confidence level (CL) of greater than 1% were retained for 
further analysis. The covariance matrix for these fits was 
optimized in a previous study [36], and the confidence- 
level distributions were checked for the present kinemat- 
ics and found to be very flat. The possible combinations 
that yield a hyperon in the strong final state are A — > p7r~ 
and E+ — > pn°. 

Figure 5 shows the invariant mass squared distribu- 
tions M 2 (p7r°) against M 2 (p7r~), summed over all kaon 
angles in one particular energy bin. For each bin in 
energy and angle, fits were done to the projections of 
M 2 (p7r°) and M 2 (p7r~) with Gaussians and a second or- 
der polynomial background. Figure 6 shows a represen- 
tative example of the fits to the A and E + peaks in a 
single energy and angle bin. After projecting and fitting 
the A and E + peaks, a region of ±3cr was chosen around 
each peak as the signal region. For further analysis of 
events with a E + , the overlap region with the A was ex- 
cluded so that there was no A distribution underneath 
the E + events. On the other hand, because the A peak is 
very narrow (approximately 1.3 MeV/c 2 across all bins 
when converted to width around the A peak, compared 
to approximately 6.3 MeV/c 2 for the E + peak), the E + 
region was not excluded from the A signal, since most 
of the A signal was within this overlap region. The re- 
maining backgrounds were removed as part of the later 
bin- by-bin yield fits. 

For these channels, a small contamination is seen in 
the projection of the invariant mass squared of the ir~ 
and 7T°, which comes from the decay K ~ — > 7r~7r°. Since 
our main channel of interest, the A(1405), is below the 
K~p threshold, we did not cut away this contamination, 
but removed it later by background subtraction. 

After these steps, we arrive at the data set of the strong 
final states of K+Ktt° and if+E+Tr". Fi gures 7 and 8 
show the invariant masses of Ytt against K + ir, where Y 
and 7t are the ground state hyperon and pion in each 
strong final state, respectively. In each of these Dalitz- 
like plots there are visible bands due to resonances in 
the Ytt system and K + tt systems. In the four ranges of 
W shown, one sees the shifting overlap of the hyperons 
E(1385), A(1405), and A(1520) versus the K*° and K*+ . 
In Section V the fits to extract the yields of each excited 
hyperon will be discussed. 



2. Event selection for E„ n and E 7r 

For these channels a final state of K + ir + ir~ with a 
missing neutron was required. A kinematic fit to the 
missing neutron mass was applied to the selected events, 
retaining those with a confidence level greater than 1%. 
Again, there are two possible hyperon combinations, 
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for a given energy bin, with projections. The bands corresponding to the A 
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FIG. 6. (Color online) Fits to the ground state hyperons (a) 
A and (b) E + for a single bin in energy and angle. The data 
of each invariant mass squared are shown as the histograms, 
and the fits are shown as the solid curve (total), dashed curve 
(Gaussian), and dot-dashed curve (background). The outer 
dotted vertical lines show the range of the fits, the inner lines 
show ±3cr around the peaks, which is used to define the signal 
events. 



— > Ti7r^~ , and these correspond to the bands shown 
in Fig. 9, where the plot of M 2 (n7r + ) versus M 2 (n7r~) is 
shown for a particular energy bin. 

To isolate the events for S + (£ _ ), we projected the dis- 
tributions onto M 2 (?i7r + ) (M 2 (n7r - )) and fit the hyperon 
peaks with a Gaussian and second order background 
polynomial. Examples are shown in Fig. 10. A region 
of ±2(7 around each peak was chosen as the signal, and 
the overlap region of the two peaks was excluded from 
each signal. Also, since there is a band corresponding to 
Kg — > 7r + 7r~ events seen in Fig. 9, we followed a simi- 
lar procedure for M 2 (7r + 7r _ ) and excluded events within 
±2cT of the Kg peak also. Figure 11 shows the invariant 
mass combinations of E~7r + versus K+TT+. Note that in 
the JC + 7r + E~ final state, the combination of K + ir + has 
no resonant structure. After this selection of events, 
the strong final states of K + T i ± tt^ z were in hand. 



3. Event selection for S°7r° 

For the £°7r° channel, the reaction is 7p — > K + T,°ir°, 
with E° — > 7A, and A — » pir~ . In this case, we were 
unable to detect the ir° as well as the 7 from the S° decay, 
therefore making a kinematic fit impossible. Instead, we 
fitted the missing mass squared (MM 2 ) with a Gaussian 
peak for the ir° and a second order polynomial for the 
signal region, and required that MM 2 be more than 3a 
above the ir° peak. Two examples of the selection of 
MM 2 are shown in Fig. 12, where the selection ranges 
are shown by the dashed vertical lines. 

To select the S° events from this channel, the invari- 




FIG. 7. (Color online) M(A7r°) versus M(K + n°) for four energy bins increasing from (a) to (d). A clear horizontal band 
corresponding to the E(1385) is seen, as well as a vertical band corresponding to the K . The contours, as well as the dashed 
lines, show the kinematic boundaries allowed in that energy bin. The blue dashed lines show the masses of each resonance 
(E(1385), K* + ) as given by the PDG, while the vertical dashed lines show where the K + ir° invariant mass is Mo ±T, and Mo 
and r are the mass and width of the K* + as given by the PDG. 



ant mass squared of the p and ir~ was plotted for each 
bin, and a fit with a Gaussian and a second order polyno- 
mial background was performed. The ±3cr region around 
the A peak was retained. The missing mass off the K + 
then gave the E°7r° line shape. For the strong final state 
K + YPir a , there are possible hyperon as well as K* + reso- 
nances, and in Section V the extraction of hyperon events 
is discussed. 

Table I shows the number of events for each channel 
after the selections shown in this section. Further selec- 
tions to isolate the states of interest will be shown later. 



IV. ACCEPTANCE AND NORMALIZATION 

To understand and correct for the CLAS detector ac- 
ceptance, a large number of Monte Carlo (MC) events 
were processed using the GEANT-based standard CLAS 
simulation package GSIM. After generating the events 
of interest, the events were passed through the detec- 
tor simulation, and the momenta were smeared to match 
the data. An earlier detailed analysis of the glla data 



showed [40] that the trigger condition for this run was not 
ideally simulated, so an ad hoc trigger efficiency correc- 
tion of ~ 5% was applied depending on the event kine- 
matics. After all corrections were made, the simulated 
events were passed through the same analysis procedures 
as the data. 

One final correction was applied for the events of inter- 
est that had a A in the strong final state. As mentioned 
in Section II, the hardware trigger for this run required 
that two particles register hits in separate sectors of the 
Start Counter. In the case of an event involving a A 
(ct = 7.89 cm [25]), there was a small probability of the 
A decaying outside of the Start Counter, and this detail 
of the trigger was not simulated in software. To rem- 
edy this, events in the simulation were removed based on 
whether the secondary vertex was geometrically outside 
of the Start Counter. The effect of this correction was 
stronger at higher energies, and for events with the kaon 
going backward in the center-of-mass frame, so that for 
most bins the correction was less than ~ 3%, while for 
some bins it was as high as 10%. For the other ground- 
state hyperons E + and £~ (cr = 2.404 cm and 4.434 cm, 
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FIG. 8. (Color online) M(£ + 7r~) versus M(K + ir~) for four bins of increasing energy from (a) to (d). Clear horizontal bands 
corresponding to the A(1405), A(1520), and higher Y* resonances are seen, as well as a vertical band corresponding to the 
K*° . The contours as well as the dashed lines show the kinematic boundaries allowed in that energy bin. The blue dashed lines 
show the masses of each resonance (A(1405), A(1520), K*) as given by the PDG, while the vertical dashed lines show where 
the K + -k~ invariant mass is Mo ± T, and Mo and T are the mass and width of the K*° as given by the PDG. 



respectively), the effect of the T, ± decaying beyond the 
Start Counter was found to be negligible. 

The photon flux in each energy bin was determined 
so that differential cross sections could be computed. 
This was done using the CLAS-standard method based 
on counting out-of-time electrons in the photon tagger 
within well-defined time windows, and correcting for the 
measured ~ 70% transmission of photons from the tag- 
ger to the physics target. Other corrections were made 
to handle tagger counters not in the primary trigger and 
to account for the measured ~ 85% data acquisition live- 
time for this data set. 



match the data. 

As mentioned at the beginning of Section III, we ex- 
tracted the S°(1385) yield in the dominant A7r° de- 
cay channel, and with the appropriate acceptance and 
branching fraction ratios scaled this down to determine 
the background yields in the IjTT^ channels. (Note that 
the E°7T channel does not result from S°(1385) decay 
because the isospin coupling coefficient vanishes.) Thus, 
the E°(1385) yield to Yiir was always known from indi- 
rect measurement within any single bin of center-of-mass 
energy and angle. For this reason, we first discuss ex- 
tracting the E°(1385) — > Air events, and then move on 
to the A(1405) yields in each T,n decay channel. 



YIELD EXTRACTION OF EXCITED 
HYPERONS 



Att u 



Our method of extracting the strong final state yields 
used simulations of the signal reaction of interest and of 
the background reactions in each channel. A fit in the 
excited hyperon spectrum was performed independently 
in each bin of center-of-mass energy and kaon angle to 



For the strong final state of -f£' + A7r , large samples of 
Monte Carlo events for the reactions 7 + p — » K + + 
E°(1385) and 7 +p — » K* + + A were generated and pro- 
cessed. For each bin in center-of-mass energy and kaon 
angle, the data and Monte Carlo events were kinemati- 
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FIG. 9. (Color online) M 2 (n7r + ) against M 2 (n7r ) for a single energy bin, with projections. The bands corresponding to the 
£~ (vertical) and E + (horizontal) are clearly seen. The faint diagonal band corresponds to the events of nKg — > mv + ir~ . 



cally fit and plotted as the missing mass from the K + , 
which is equivalent to the invariant mass of the A and 
7T°. A fit to the data with these Monte Carlo templates 
was performed by scaling each Monte Carlo template by 
an overall factor. Figure 13 shows the fit results for some 
of these bins. In all bins the K* + A channel contributes 
as a smooth background. 

A peculiarity was noticed that the £(1385) line shape 
could not be fit well using a relativistic Breit-Wigner 
function with a mass-dependent width. Rather, a non- 
relativistic Breit-Wigner function with width indepen- 
dent of mass was seen to match the data much better. 
The line shapes used as input to the CLAS Monte Carlo 
for counts dC x (to) / dm as a function of Ytt mass m were 

dC aon _ Iel (m) r /27r 

dm ~ (too - m) 2 + (r /2) 2 1 ' 

and 

dCrejjm) ^ (2/ir)mm r(q) 
dm (toq — m 2 ) 2 + (m,QT(q)) 2 

for the non- relativistic and relativistic cases, respectively. 
The mass-dependent width was T(q) = F a (q/q ) 2L+1 , in 
which q (go) is the breakup momentum of the A7r° or 
£-7r system in the resonance rest frame at mass m (mo). 
The orbital angular momentum in this case is L = 1. 
Figure 14 shows a comparison of the fit results of the 
£°(1385) peak using Monte Carlo templates generated 
with the forms of Eqs. 5 and 6. Clearly the relativistic 
Breit-Wigner template is not able to fit the data well, 
essentially because the g 3 factor in the numerator sup- 
presses the yield near threshold too much. Therefore 



for our present purpose, we used the very simple non- 
relativistic Breit-Wigner form for fitting the £°(1385) 
data in each bin. 

The reason why the non-relativistic Breit-Wigner form 
fits better to the data is not clear, but we note that previ- 
ous experiments that determined the £(1385) mass and 
width based on hadronic reactions [25], also had diffi- 
culties in fitting to a relativistic P-wave Breit-Wigner 
line shape, and tested non-relativistic forms with mass- 
independent widths [41-43]. Because these papers mea- 
sured the charged £(1385) line shapes, where leakage due 
to the A(1405) or other A* states cannot occur, this seems 
to be an inherent feature of the £(1385), and not due to 
some unaccounted leakage in our data. Furthermore, the 
effect is seen across all of our energy bins, even when 
below the nominal K* + threshold or when kinematically 
separated from the K* + . Therefore, we conclude that 
this effect is not due to interference with the K* + 1 . 

After the yields of the £°(1385) were extracted in 
each bin of center-of-mass energy and angle, the differen- 
tial cross sections were calculated using the acceptance 
based on simulations and the photon flux normalization. 
The £°(1385) differential cross section results will be dis- 
cussed in a separate paper [32], along with those for the 
A(1405) and A(1520). In this paper we focus on extract- 
ing the yields for the A(1405), for which the yields of the 
£°(1385) decaying to £ ± 7r =F are necessary. 



1 Refs. [41-43] used the reaction K + p — > A + 7r+ + n , where 
there is no K* + background. 
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FIG. 10. (Color online) Fits to the ground state hyperons (a) 
E + and (b) E~ for a single bin in energy and angle. The data 
of each invariant mass squared are shown as the histograms, 
and the fits are shown as the solid curve (total), dashed curve 
(Gaussian), and dot-dashed curve (background). The outer 
dotted vertical lines show the range of the fit, and the inner 
lines show ±2cr around the peaks used to select events. The 
opposing signal region was excluded, as well as the ±2cr peak 
around the Kg. 



For each of the charged Y,ir channels, the acceptance- 
corrected yield of the An channel (BR = 87.0%) was 
scaled down by using the branching ratio (BR — 11.7%) 
and acceptance for each bin. Because the E°(1385) yield 
was based on a measurement of the Air channel, it was 
not allowed to vary when extracting the yields of the 
A(1405). 



B. E + 7T" 

We next focus on the £+7r~ channel, although the 
other E7T channels are quite similar in procedure. As can 
be seen from the plots of M(E + 7r~) versus M(K + tt~) in 
Fig. 8, there are contributions from the A(1405), A(1520), 
and other excited hyperon states, as well as from the 
K*°. We model each of these contributions separately 
with Monte Carlo event templates. Each template is 
generated according to a relativistic Breit-Wigncr form 
with its resonance mass Mq and width Tq taken from 
the PDG [25]. We assumed a mass-dependent width of 
T(M) = T (q/q ) 2L+1 , where q (q ) is the breakup mo- 
mentum of the daughter particles in the resonance rest 
frame at mass M (Mq) with L the orbital angular mo- 



mentum. In the fitting procedure, only the normalization 
of each template was allowed to change to get the best 
agreement with the data. For the £°(1385) contribution, 
the yield was fixed by the A7r° channel discussed above, 
and therefore the yield was not allowed to vary. 

Figure 15 shows a fit result for the £ + 7r - invariant 
mass spectrum using the above templates for a single bin 
in center-of-mass energy and angle, along with a back- 
ground Breit-Wigner function that fits the Y* resonance 
around 1670 MeV/c 2 . Since our goal is to extract the 
A(1405) line shape in the most model-independent way, 
we start with a relativistic Breit-Wigner form based on 
the PDG values of mass and width for the A(1405), and 
this is shown in Fig. 15(a) as the red points. The fit 
is inadequate around the A(1405) region, showing that 
a simple Breit-Wigner function is not able to describe 
the data well. For this reason, the template form of the 
A(1405) was modified in an iterative way, as explained 
below. 

Once an initial fit was obtained, we subtracted inco- 
herently the contributions due to the E°(1385), A(1520), 
K*°, and F*(1670) so that the only remaining contri- 
bution was from what should be the A(1405). We call 
this the residual distribution for the A(1405). Because 
this residual distribution is the best measure for the raw 
A(1405) yield, we applied an acceptance correction based 
on the Monte Carlo simulation of CLAS. A large number 
of events was generated flat in the K + Yj + ii~ three-body 
phase space, and the residual distribution was corrected 
as a function of the £ + 7r~ invariant mass based on the 
acceptance of these events. After acceptance correction, 
the true line shape of the A(1405) was obtained for each 
energy and angle bin. 

As noted above, the A(1405) was not adequately 
described by the initial template, so we used the 
acceptance-corrected line shape obtained with the above 
procedure to iterate the Monte Carlo template for the 
A(1405) region. The iteration process made use of data 
summed over all kaon angles within each energy bin. Fig- 
ure 15(b) shows the third and final iteration. Note that 
the total fit is now closer to the data, and we see how 
the iteration converged to stable line shapes based on 
the data. Since the residual is determined by subtracting 
off components such as the £(1385) and the A(1520), the 
residual shapes do not depend strongly on the exact tem- 
plate shape we used for the A(1405). The residual that 
was obtained from the fit using this final template was 
acceptance-corrected and normalized to the photon flux, 
yielding our intermediate result for d 2 a/dmd cos in 
bins of energy, angle, and £ + 7r~~ mass m. 

The procedure for the £+7r~ channel was exactly the 
same as for the S+7r~ channel, because the physics is 
identical except for the final decay of E + — > mr + . Line 
shapes were obtained in each energy and kaon angle bin. 
By comparing the two E + 7r~ channels we were able to 
check our results, as will be shown in Section VI. 



12 




FIG. 11. (Color online) M(E~7r + ) versus M(K + tt + ) for four energy bins increasing from (a) to (d). Clear horizontal bands 
corresponding to the A(1405), A(1520), and higher Y* resonances are seen. The blue dashed lines show the nominal masses 
of the A(1405) and A(1520) from the PDG. Note that in this channel the combination of K + ir + shows no resonant structure. 
The contours as well as the dashed lines show the kinematic boundaries allowed in that energy bin. 



C. E-7T + 

For the E~7r + channel, we followed the same proce- 
dures as above, but in this case, with the strong final 
state of K + Y,~tt + , the Kir combination is exotic, and 
therefore we expect no resonance. However, to accom- 
modate the broadly distributed events seen in Fig. 11, a 
phase-space distribution of K + E _ 7r + was generated, and 
this was used as a fit component. The line shapes were 
iterated as before and then acceptance-corrected. 

An interesting feature of this channel is the presence of 
the y*(1670), which shows up much more strongly com- 
pared to the other E-7T channels, as seen in Fig. 11. The 
PDG lists several candidate resonances in this region, but 
we have not made an effort to further identify this state. 

D. EV 

For the remaining E°7r° channel, we did fits to the K + 
missing mass distribution similar to the previous cases, 
but since the E°(1385) cannot decay to E°7r due to the 
vanishing isospin factor, there is no E°(1385) contribu- 



tion. The fits were performed with templates for A(1405), 
A(1520), and K*+YP . As the F*(1670) region does not 
show any prominent peaks, the Breit-Wigner function for 
y*(1670) was not used. 

In summary, all E-7T channels were isolated in order to 
extract the line shape of the A(1405) region based on 
fits to the data. The line shape templates for the region 
of interest were generated in an iterative way using the 
data for each channel independently, and in all cases the 
results showed convergence after several iterations. 



VI. LINE SHAPE RESULTS 

The E7T mass distributions or line shapes, 
da I 'd cos C ^+ dm, were obtained in each bin of center- 
of-mass energy and kaon production angle, but due to 
limited statistics we have summed over all angles within 
each energy bin to obtain a single line shape, da /dm, for 
each energy bin. Alternatively, we can sum over mass 
to obtain da/dcos9 c K f- , the differential cross section. 
These results will be shown in a separate paper. Here 
we compare the results of the two E + 7r - channels for 



13 



400 
300 



c 

8 200 



a 

3 

o 



100 


800 



600 



400 



200 





1:1 1 

2.05<W<2|.15 (GeV) 


F/ 1 

7 1 

— jd / \V 


O.7O<cos0^f<O.8O 

— 


0.02 


0.04 „ 0.06 „„ 0.08 0.1 
MM 2 (GeV/c 2 ) 2 




1:1 1 

2.35<W<2.45 (GeV) 


-A 

-/ 1 


0.70<cos9^<0.80 _ 

— 


j. \ 




0.02 


0.04 0.06 „„ 0.08 0.1 
MM 2 (GeV/c 2 ) 2 



FIG. 12. (Color online) Examples of the 7p — > _K" + p7r~(V) 
missing mass squared (MM 2 ) spectrum for selected kinematic 
bins. The vertical dashed lines show the selection range for 
the S°7r° channel between the 7r° and 2tv° limits. 



consistency, then proceed to a comparison of all three 
E7T channels. 



A. Line shape results for S + 7r channels 

We begin by studying the two channels S+7r~ and 
which share the same strong final state, and dif- 
fer only in the decay of the S + . Comparing these two 
channels gives a measure of the reconstruction accuracy 
of the analysis. Figure 16 shows a comparison of the 
line shapes obtained for each £ + 7r~. channel. The in- 
ner error bars are the combined statistical uncertainty 
of the data and of the Monte Carlo samples that were 
used in the background subtraction. Our fits to the raw 
invariant mass spectra using Monte Carlo templates did 
not always perfectly reproduce the data, even after iter- 
ating. To account for this possible systematic error in 
our analysis, we summed the data within each energy 
bin over all kaon angles and compared to the summed 
fit result. Any discrepancy in each mass bin was taken 
as an additional uncertainty, and a portion was added 
in quadrature with the statistical errors above. These 
are shown as the outer error bars in Fig. 16. Thus, the 
outer errors bars represent the combined point-to-point 
statistical and systematic uncertainty. 

Beyond these estimated uncertainties on single decay 
modes, any large, possibly-nonstatistical difference be- 
tween the two measured £ + 7r~ modes could also signal 
a systematic discrepancy in the analysis. Therefore, for 
each mass bin, we took the difference of the two measured 



values and subtracted the summed errors in quadrature, 
obtaining a mass-dependent error that estimates this sys- 
tematic discrepancy. The shaded histogram at the bot- 
tom of the plots shows these uncertainties when the dif- 
ference of the two measured points is larger than the 
sum of the two errors. The agreement between the two 
decay mode reconstruction channels is generally good. 
The average of these two measurements will be used in 
the subsequent comparisons with the other charge decay 
modes. 

In all cases the £ + 7r~ mass distribution clearly peaks 
at a mass of around 1420 MeV/c 2 which is higher than 
the nominal mass of the A(1405) at 1405.1 MeV/c 2 listed 
by the PDG [25] . We also note the sharp drop or break of 
the mass distributions at the NK threshold near 1.435 
GeV/c 2 , which is a signature of the opening of a new 
threshold for S'-wave resonances. This will be discussed 
in Section IX. 



B. Line shape results for all Y>n channels 

Our main results [44], the line shape comparison for 
all three T,ir channels, is shown in Fig. 17. As noted, the 
£ + 7r~ channel is the weighted average of the two mea- 
sured final states. The E°7r° channel and £~7r + channels 
are again shown with inner and outer error bars, where 
the inner bars are statistical, and the outer bars include 
the estimated residual discrepancy in the fits added in 
quadrature to the inner bars. For each of nine bins 
in invariant energy W, we show the T,tt mass distribu- 
tion in each of three charge states. The data have been 
summed over the full range of measured kaon production 
angles. The large-angle cut-offs were not quite identi- 
cal for all charge states because of differing acceptances, 
but since the cross sections get very small at large angles 
(cos 6f£+' < —0.5) we can neglect these differences. 

For all energies, it is evident that the line shapes dif- 
fer markedly between charge states; in some regions they 
differ by well over 5a. This occurs far away from the in- 
dicated reaction thresholds, making it unlikely that the 
effects are due to mere mass differences. None of the 
mass distributions are reproduced by the simple relativis- 
tic Breit-Wigner line shape with PDG-given centroid and 
width. The £ + 7r~ channel peaks at a higher mass than 
the £ - 7r + channel, while having a width that is signifi- 
cantly smaller. The charge-dependence of the mass dis- 
tributions is largest for W between 2.0 and 2.4 GeV. For 
W approaching 2.8 GeV the mass distributions tend to 
merge together. This hints that whatever / ^ coherent 
admixture of isospin states is at work here, it fades away 
at higher total energy. Our own fit to the line shapes 
to extract our best estimates for the mass and width 
of the A(1405) and other structures causing this charge- 
dependence of the mass distributions will be shown in 
Section IX. 

Comparing our line shape results to the prediction of 
Nacher et al. [7] computed in a chiral unitary model ap- 
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FIG. 13. (Color online) Sample fit results of the strong final state of K + Air" . The events are plotted versus the missing mass 
from the K + , which is equivalent to the invariant mass of the An system. The data are shown by the black crosses, while the 
E°(1385) signal Monte Carlo and the K* + A background are shown by the red crosses and green circles, respectively. The sum 
of the Monte Carlo templates are shown by the solid magenta line. 




FIG. 14. (Color online) Sample invariant mass spectra for W = 2.6 GeV and cos (9^+' = .65 showing the E°(1385) peak 
using a Monte Carlo template based on (a) relativistic Breit-Wigner (mass-dependent width) (b) non-relativistic Breit-Wigner 
(mass-independent width). The fit with the relativistic Breit-Wigner form clearly does not fit the data well. 




FIG. 15. (Color online) Fit result to the strong final state of E + 7r~ before and after Monte Carlo iteration, for W = 2.5 GeV, 
cosd c K ™' = .35, as a function of the E + 7r~ invariant mass. The data are shown with black crosses, (a) Before iteration. Each 
Monte Carlo template and the Breit-Wigner function for the Y*(1670) is shown by a separate color. The total simulation is 
shown in cyan, (b) After iterations of the A(1405) template. 



15 



> 

o 



X5 




Z + 71" — > p JC° JI" 
Z + 71" — > tl 7I + 71" 

weighted average 
PDG Breit-Wigner 



1.5 



1 7t Invariant Mass (GeV/c ) 

(a) 



> 

o 


-a 



1.5 
1 

0.5 

0, 



2.25<W<2.35 (GeV) /"\ 



E + 7T" threshold 



'1.3 



E + % — > p 71° 7T." 
£ + 7T." — > 11 7C + 71" 

weighted average 
PDG Breit-Wigner 




E 7t Invariant Mass (GeV/c ) 

(b) 



> 

o 



T3 



0.5 



2.45<W<2.55 (GeV) 



£ + 7t" threshold 



u 1.3 



Z + 7C" — > p 71° 71" 
Z + 7C" -» tl 7I + 71" 

weighted average 
PDG Breit-Wigner 



/ Wl-k i* - 1 i 

fit * iFr;-,-,.i. s . ! 

'm i — r 

1.4 1.5 
E Jt Invariant Mass (GeV/c 2 ) 

(c) 



0.6 



> 

u 

O 0.4 



s 

T3 

D 
-a 



0.2 



2.65<W<2.75 (GeV) / 



L + 7T" threshold 



u 1.3 



iitr 



1.4 



E + 71" — > p 71° 71" 
£ + 71" ^ n 7X + 7t" 

weighted average 
PDG Breit-Wigner 



1.5 



E 7t Invariant Mass (GeV/c ) 

(d) 



FIG. 16. (Color online) Line shape results for the two E + 7T~ channels. The Ep tt - channel is shown with light magenta open 
circles, while the E^7T~ channel is shown with light blue triangles. The weighted average of the two line shapes is taken as the 
final E + 7r _ line shape, and is shown as red downward triangles. The dashed line represents a relativistic Breit-Wigner function 
with a mass-dependent width, with the mass and width taken from the PDG. The blue-hatched histogram at the bottom shows 
the averaged estimated systematic discrepancy between the two reconstructed decay modes. 



proach, we see in Fig. 18 that they are indeed different 
for each Utt channel. In the chiral unitary theory this 
was explained as an / = 1 amplitude interfering with the 
1 = A(1405) amplitude in such a way that the S + 7r _ 
and E _ 7T + channels were shifted in opposite directions 
due to the interference term. The model curves were 
computed for E 7 = 1.7 GeV, but we compare with our 
results at E~ = 1.88 GeV since our statistics are better 
there. The model calculation uses a Weinberg- Tomozawa 
contact interaction that is energy and angle independent, 
allowing us to compare the model to the data in any en- 
ergy bin. In our results it is the E + 7r~ channel that 
is shifted to higher mass with a narrower width, and 
the E~7T + channel is smaller and wider, in contrast to 
the model calculation. Also, the model curves have been 
scaled down by a factor of 2.0 to match the data, sug- 
gesting that the model overestimates the strength of the 
photocouplings by that amount. In Section IX we will 
make our own phenomenological isospin decomposition 
to find a plausible explanation of what is seen. 

The other existing prediction for the mass distribution 
of the Yj-k final states is that of Lutz and Soyeur [11]. 
In their so-called double kaon pole model, the combined 
effects of the E(1385) and the A(1405) were considered, 
and this produced some variation among the three charge 



combinations we have presented. However, as has been 
discussed, we subtracted off the effect of the E(1385) and 
still are left with a substantial variation in the three fi- 
nal states. We do not compare our results directly to 
theirs because they are qualitatively similar in shape to 
those of Ref.[7], and also because they are about a fac- 
tor of four too large in cross section, indicating a serious 
quantitative discrepancy when comparing to our results. 



VII. SYSTEMATIC UNCERTAINTIES AND 
TESTS 

A. Overall systematics of the run 

For systematic uncertainties, there were global contri- 
butions from the yield extraction, acceptance corrections, 
flux normalization, and the line shape fitting procedure. 
The main cuts that influenced the yield extraction were 
the ATOF cuts, the confidence level cuts in the kinematic 
fit, and the selection of intermediate the ground-state hy- 
peron. All of these cuts were varied within each bin of 
center-of-mass energy and angle, and the total yields were 
checked for any differences due to the cuts. Variation 
in the ATOF width by 0.2 ns changed the acceptance- 
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FIG. 18. (Color online) Mass distributions at W = 2.10 GeV 
and _E 7 = 1.88 GeV in comparison to the model of Nacher et 
al. [7] scaled down by a factor of 2.0. The E + 7r~ channel is 
shown as red circles and the red dot-dashed line; the E°7r° 
channel is shown as the blue squares and the blue dashed 
line; the E~7r + channel is shown as the green triangles and 
the green solid line. The dashed vertical colored lines at the 
left side show the reaction thresholds, and the vertical dashed 
lines at 1.405 GeV and 1.437 GeV mark the nominal centroid 
and the NK thresholds, respectively. The error bars on the 
data points are combined statistical and point-to-point sys- 
tematic uncertainties. 



corrected yield between 2-6% in each bin. Changing the 
confidence level from the nominal 1% to 10%, changed 
the acceptance-corrected yields by 3-12% in each bin. 
For the majority of the bins, the final yields changed by 
less than 4% for the ATOF cuts, less than 7% for the con- 
fidence level cuts, and less than 2% for the ground-state 
hyperon selection. 

Stability of the normalization was monitored through- 
out the run. The fluctuations in target density were de- 
termined to be a negligible 0.11%, while the target length 
was measured to 0.125%. The photon normalization was 
examined on an hour-by-hour basis by measuring the w 
production yields [37], and the uncertainty for the nor- 
malization was determined to be 7.3%. The live-time 
correction that was necessary to determine the photon 
flux introduced an additional uncertainty of 3%, and the 
photon transmission efficiency added 0.5%, so that the 
total uncertainty for the photon normalization was 7.9. 

For the final systematic uncertainty all of the above 
global uncertainties were added in quadrature to yield a 
final value of 11.6%. A summary of each uncertainty is 
shown in Table II. 

The mass resolution of the line shape results was inves- 
tigated by generating Monte Carlo samples of zero width 
centered at E-7T invariant masses of 1.406, 1.450, and 1.500 
GeV/c 2 . Because the E-7T invariant mass is equivalent to 
the missing mass off the detected K + , the E-7T mass res- 
olution was related to the momentum resolution of the 
K + . However, kinematic fitting of most of the channels 



TABLE II. The global systematic uncertainties in the exper- 
iment. They arise from yield extraction, acceptance calcu- 
lation, target characteristics, photon flux normalization, and 
branching ratios [25]. The total was calculated by summing 
all in quadrature. 



Source 


Value 


ATOF cuts 


2-6% 


Confidence level on kinematic fit 


3-12% 


Selection of intermediate hyperons 


2-3% 


Target density 


0.11% 


Target length 


0.125% 


Photon normalization 


7.3% 


Live-time correction 


3% 


Photon transmission efficiency 


0.5% 


£(1385) -> E7r, Att 


1.5% 


A — > pn~ 


0.5% 


E+ -> pn°,nn + 


0.30% 


E~ — > nn~ 


0.005% 


Total 


11.6% 



improved the overall mass resolution. For all generated 
events in all bins of center-of-mass energy and angle, the 
E-7T invariant mass for the accepted events was fit with a 
Gaussian to determine the resolution. This showed that 
for the lower energy W bins, the resolution (a) was better 
than 6 MeV/c 2 , while for the higher energy bins it was 
up to 8 MeV/c 2 , with worse resolution in the backward 
kaon angles, where the CLAS magnetic field is weaker. 
Without the kinematic fit, as in the E°7r° results, the 
mass resolution averaged about 2 MeV wider at high W 
and large angles. No shift of the center of the Gaussian 
larger than 1 MeV/c 2 was seen. Since our results are 
shown with 5 MeV/c 2 bins, the mass resolution of the 
line shapes is one to two bins. We also remark that the 
absolute mass accuracy of the experiment for hyperons 
such as the E(1385) and the A(1520), and of meson states 
in this mass range is < 1 MeV/c 2 . 

B. Removal of K* 

A concern in the photoproduction line shape analysis 
of the A(1405) region is the effect that the K* may have. 
As seen in Fig. 8, the A(1405) has a kinematic overlap 
with the K* in the strong final-state phase space, so that 
the difference in line shapes seen in the various E7r chan- 
nels could be due to interference with the K* . Below we 
argue that this is not the case. 

We measure the line shape of the A(1405) in bins of 
center-of-mass energy, W, and the kinematic overlap of 
the K* depends strongly on this energy. Fig. 8 shows 
no K* overlap at low W below the K* threshold, strong 
overlap at intermediate W and again no overlap at high 
W. 

We tested for the presence of K* interference by cut- 
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ting out regions of K + ir~ invariant mass centered around 
the K* mass and in multiples of T/2, where T is the width 
of the K* listed in the PDG. Figure 19(a) shows the ef- 
fect of each K* rejection cut up to ±|T at W — 2.0 GeV. 
For each cut, we reprocessed all of the Monte Carlo sam- 
ples of the other channels used in the template fit, redid 
the fit, and applied acceptance corrections. If there were 
any interference between the K* and the A(1405), we 
would expect it to be strongest in the region where the 
K* is strongest, whereas our results show that even with 
the overlap region removed, the final result is remark- 
ably unchanged by this drastic removal. In Fig. 19(b), 
the line shape changes significantly only for the cut at 
±|r (green downward triangles), but this is simply due 
to the loss of phase space and acceptance, since the cut 
removes about 300 MeV/c 2 of K + ir~ invariant mass cen- 
tered around the K* . This is reflected in the Dalitz-like 
plot of the strong final state in Fig. 8, where the bound- 
aries of ±r around the K* mass are shown as vertical 
dashed lines. Figures 19(c) and 19(d) are more evidence 
of the insensitivity of the S7r mass distributions to the 
K* . A similar study was done to test for possible co- 
herent interference between K* + A and if + E(1385), and 
again no such effect was detected. 

We conclude that, although we cannot completely rule 
out interference due to the K* , our results are not signif- 
icantly altered even when we apply a drastic cut on the 
K* region, thereby removing most of its strength. Be- 
cause the photoproduction line shape of the A(1405) is 
not known to any accuracy, we do not attempt any fur- 
ther analysis of the interference with the K*. We antic- 
ipate our measurement will further stimulate theoretical 
interest in this state, and with more theoretical input, a 
more elaborate analysis may be possible in future exper- 
iments. 



that we will consider in order to gain a reasonable rep- 
resentation of the experimental results. The first arises 
from the channel-coupling between the detected T>ir fi- 
nal state and the undetected NK final state. This is 
done by using a Flatte-like formalism [45] to enforce two- 
channel unitarity and analyticity of the production am- 
plitude. The second arises because we find that the dif- 
ferent charge states have markedly different mass distri- 
butions, implying that amplitudes other than 1 = must 
participate in the reaction mechanism. This is treated by 
including coherent 1=1 amplitudes that interfere with 
the 1 = amplitude. 

Since the electromagnetic interaction does not conserve 
isospin, the initial 72? state in this reaction can have both 
I = 1/2 or 3/2 character, and will lead to a final K + (Ett) 
state wherein the £71 system is in a superposition of 
It,tt = 0,1, and 2 states. The three measured Y,ir fi- 
nal states all have their third component of isospin, 
equal to zero. If we denote the isospin state of the system 
as I/ttEi^s)' we can wr ite each of the three measured 
final charge combinations using Clebsch-Gordon coeffi- 
cients as 
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(7) 

(8) 
(9) 



Let tj be the complex matrix element that takes the 
initial jp state via a transition operator to the final 
state that contains the kaon and the T,tt system in the 
/stt = I state, so that 



VIII. MODEL FOR ISOSPIN DECOMPOSITION 

We have shown that the line shapes of the E7r final 
states are far from those of a simple Breit-Wigner form. 
Indeed, there are two main modifications to the picture 
of a simple single resonance for the A(1405) mass region 



\t^ = \(I,0\fW\jp}\' 



(10) 



The magnitude-squared matrix element for creating 
a particular charged final-state pair, T^a-^b (a, b € 
{H — , 00, — h}), can then be obtained by combining these 
expressions. For example, the probability of populating 
the |7r~£ + ) state is proportional to 



KTr-E+lfW+fW+f^lTp)! 5 



;N 5 



V6 



|*o||*i| cos< 



'01 



1*1 1 1*2 1 COS</> 



V2 



12 



1*0 1 1*2 1 COS 002, 



(11) 

(12) 



in which the real relative phases between the three isospin 
amplitudes are </>oi(wi), fyxzijn), and 0o2( m )- The other 
two charge combinations have similar forms. Thus there 
are five real parameters, assuming one phase is set to 
zero. We expect the matrix element to have the kine- 
matic dependence = T„-- s +(W,m), where W is 



the available overall center-of-mass invariant energy and 
m is the T,tt invariant mass. At a given value of W and 
m we have three measured cross sections that are pro- 
portional to the three quantities T^a^h, so we cannot de- 
termine all five numbers uniquely. 

Before going on, we chose at this point to apply the as- 
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FIG. 19. (Color online) Final results for the line shape in the E^7r channel when the K* is removed successively in steps 
of ijT, where Y is the width of the K* quoted in the PDG. The black circles represent our final results without a cut on the 
K* , while the red squares, blue triangles, and green downward triangles represent cuts of ±r/2, ±F, ±|F centered around the 
K* , respectively. 



sumption that the 1 — 2 amplitude is negligible, and that 
all of the interference in this reaction is between 1 = 
and 1 = 1 amplitudes only. This assumption is consis- 
tent with all previous work on this subject, for example 
Refs. [7, 9], in which the dynamics of the A(1405) is pre- 
sumed to be all within 1 = and/or both 1 = and 1. 
With this assumption, we can write the expressions for 
the production strength of the three E7r channels as 



The 1=1 amplitude's magnitude squared is given by 



ti 



1 21 



■s+l 



121 



I 2 — 21X^020 | 2 



(17) 



which implies that the average of the charged final states 
should be greater than or equal to the neutral final state, 
depending on the size of The interference between 
the isospin states is accessed using 



l^- s+ | 2 = ^ | 2 + \\ti\ 2 - -jg|to||*i| cos^oi, (13) 



(14) 



|2>£°| 2 = cil^l 2 ' 

l^ +s -| 2 - ^\t \ 2 + ^\h\ 2 + ^=MM cos^oi- (15) 

These relationships can be combined to show several 
things. First, the sum of the measured line shapes gives 
the sum of the 1 = and 1=1 amplitudes' squared 
magnitudes: 

ir.-E+l 2 + |2> s0 | 2 + |i; +s -| 2 = \t \ 2 + M 2 , (i6) 

that is, the interference terms cancel and we see the inco- 
herent sum of the isospin channels. The 1 = amplitude 
is proportional to the S°7T° channel alone, as per Eq. (14). 



IT; 



| 2 -|T.- S+ | 2 



1*0 1 1 *x I COS 001. 



(18) 



This equation shows how any difference between the 
charged decay modes is directly related to the interfer- 
ence of the two isospin channels. Note that </>oi is the 
mass-dependent phase between to(m) and t\{m). Apart 
from that mass dependence, we allow an arbitrary strong 
production phase for each of the amplitudes, called A</>/ 
below. 

For the production reaction 7 + p — > K + + (EV) we 
write the contribution from an amplitude of isospin I at 
fixed 7p center-of-mass energy W and £77 mass m as 



t x (m) = C/(PF)e lA ^ / B / (m), 



(19) 



where Cj{W) is a real number representing the effective 
strength of the excitation and A</5j is a corresponding 
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production phase. The Breit-Wigner amplitude has the 
form 



Bi{m) = 



Vraj?rar°(q/g K ) 2£ 
m R -m 2 - im R T tot (q) 



(20) 



where m R is the centroid of the resonance distribution, 
in this case the E7r invariant mass, is the fixed decay 
width to a given final state, and r tot (q) is the total width 
to all final states. The available momentum in the decay- 
ing hyperon center-of-mass system is called q, and in this 
frame qn is the available decay momentum at m = tur. 
In this way of writing the amplitude, the numerator has 
no phase space factor, but this will be included below 
when we write the final expression for the line shape. 

We assume that the line shape for each isospin contri- 
bution to the intermediate hyperon state is described by a 
relativistic Breit-Wigner distribution with suitable phase 
space factors and normalization. The total width of the 
resonance, T tot (q), is the sum of partial decay widths, 
but for a single decay channel designated by a "1", let 
it be the partial decay width Tj i(qi{m)). The width of 
the resonance going into a single decay mode "1" is, in 
the relativistic formulation, dependent on the mass and 
is written as 



' m 



qi(m) 
qn 



2L+1 



(21) 



where 1 denotes a fixed decay width that will be deter- 
mined by the fit, and q\{m) is the available momentum 
in this decay mode at mass m. This expression accounts 
for the increasing phase space available for the two-body 
decay across the resonance, and it forces the width to 
zero at threshold. (Later we will analytically continue q 
below threshold, however.) We will consider only L = 
or S-wave decays, as required for the odd-parity A(1405) 
decaying to a pseudo-scalar meson and an octet baryon. 

The overall coupling strength of the resonance repre- 
sented by Eq. (20) for the reaction 7+p ->■ AT++A(1405) 
is given by the parameters in Eq. (19). We take these to 
be fixed (at a given value of W) over the whole range of 
the mass distribution m. 

For several of the fits to the data (discussed below) 
we used either two I = or two 1 = 1 Breit-Wigner 
amplitudes. In all cases these amplitudes were added 
coherently. We selected A0o for the 'first' or 'dominant' 
1 = amplitude to be zero, so the other strong phases 
were determined relative to it. 

For hadronic reactions we must also consider the dy- 
namical consequences of the opening of thresholds to de- 
cay channels other than the single channel denoted "1". 
In the present situation there is the NK channel that 
opens at m t hresh = m K + m N — 1434 MeV/c 2 , which is 
within the range of the mass distribution of the E-7T sys- 
tem under study. This can significantly impact the line 
shape of the resonance. To preserve unitarity and the 
analytic form of the decay amplitude as a mass thresh- 
old is crossed, we modify the amplitude of Eq. (20) in 



a specific way. If we denote the second decay mode as 
channel "2" . then the total width of the resonance is 



"tot(m-) = T I s(q 1 (m)) + r I>2 (q 2 (in)), 



(22) 



where the second decay channel is described by the width 

2L+1 



r <n\ ro mR ( 92 ( TO ) 
1 1,2(1) = 1 1,2 



m \ qn 



(23) 



Here, q2{m) is the decay momentum available for decay 
mode "2" at mass m, and 2 is the constant factor 
for the width of this partial decay mode. Below thresh- 
old TOthreshj the momentum q2(jn) is nominally zero. 
However, in the Flatte method [45] we analytically con- 
tinue the momentum to imaginary values, denoting it as 
qj = —iq2 for m < mthresh- Furthermore, we introduce a 
Flatte parameter for the branching fraction of the decay 
modes as 

7 = r?, 2 /r? ;1 . (24) 

Below threshold for decay mode 2, the total decay width 



r tot (m) = r 7 ,i( ? i(m)) + z 7 r/,i(g 2 'M), (25) 
while above the threshold the total decay width is 



r t ot(m) = Y\{q\{m)) +jTi(q 2 (m)), 



(26) 



and these two expressions are used, respectively, in 
Eq (20). 

Apart from the overall strength Cj and phase A</>/, 
there are two free parameters in these expressions for a 
single resonance: the intrinsic width Tj 1 , and the relative 
branching fraction between decay modes 7. The fits were 
made over the whole range of energy W (in nine bins from 
2.0 to 2.8 GeV), and these two parameters were fixed to 
the same value for all W. 

The experimental results for kaon-angle integrated 
mass distributions are in the form of differential cross 
sections da a b/dm with ab € {H — ,00, — h}; the expres- 
sion for this cross section includes relevant flux and phase 
space factors. Fig. 1 illustrates how this reaction requires 
the use of three-body phase space. To arrive at it we fac- 
torize this phase space into two two-body pieces using 
standard methods [46], the first for the K + Y* hyperon 
intermediate state of mass m, and the second for the de- 
cay of this state into £71". 

The fully differential form of the cross section is 



d<Tab 



= (M 2 PK+q 

dfl^+dfl^dm (47r) 5 p lp s 



T 

-*- TT 



(27) 



where the momentum of the kaon in the overall center-of- 
mass system is Vk+ 1 the momentum in the £7r final state 
is the aforementioned q, and i/s = W. The invariant 
production amplitude is defined by Eqs. (13), (14) 
and (15). 
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In the experiment we measure the decay distribution 
of Y* —> T.TT over the full solid angle f^, so the data are 
automatically integrated over this variable. Formally, we 
take Ttts to be independent of this decay angle. The 
reaction is not "flat" with respect to kaon angle, as we 
know from measurement of the differential cross section 
da/dfl K + . But for studying the line shapes we are forced 
to integrate over kaon angle to gain enough statistics for 
the analysis. We therefore take T^s to be the kaon-angle 
averaged matrix element and integrate over £Ik+ ■ Fig- 
ure 1 shows that there is a vertex involving the photon, 
and the strength at this vertex must be proportional to 
y^a, where a is the fine structure constant. Factoring 
this out of the matrix element means the previously- 
defined fit parameters Cj become an effective strong cou- 
pling with units of v GeV. The final expression for the 
differential-in-mass cross sections is then 



da ah (he) 2 a p K +q 



dm 64-7T 3 p lp W 2 



T. 



(28) 



There is an interplay among the phase space factors in 
front of the matrix element. For a given invariant energy 
W, p lp is determined. But the possible ranges of p K + 
and m are also limited, so the larger px+ becomes, the 
smaller m, and therefore the smaller q must be. 

In addition to the coherent sum of the isospin com- 
ponents of the line shapes, it was necessary to include a 
linear background function under each of the E-7T mass 
distributions. This sloping background was introduced 
to represent less-than-perfect subtraction of the back- 
grounds due, for example, to K* production or tails of 
higher mass hyperons. The need for such a background 
parametrization is seen in the data, which show that in 
several mass distributions the trend at the high-mass end 
of the scales is not toward zero, but rather to a constant 
or even a rising slope. The problem was mainly with 
the E°7T final state, the one for which it was not pos- 
sible to make a direct experimental measurement of the 
K* background, and we had to rely on the Monte Carlo 
model alone. The slopes of the backgrounds were not fit 
parameters, but were matched to the differential mass 
distributions at 1.6 GeV. 



IX. ISOSPIN DECOMPOSITION 

We can now take the mass distributions found in this 
analysis and separate the information from the three 
charge combinations in the E7r final states according to 
1 = and 1=1 components. This is crucial toward 
the goal of understanding the contribution from the true 
A(1405), which is by definition 1 = 0, and anything else 
happening in the reaction mechanism. 

We found that fitting two 1 = amplitudes to just 
the S°7r° data led to a very good fit after including the 
Flatte channel coupling [47] . A "two-pole" explanation of 
the A(1405) would favor such a result. The centroids and 
widths of the 1 = states remained stable when an I = 1 



amplitude was added to include the £ + 7r~ and £~7r + 
final state combinations. However, it was found that a 
much better fit could be obtained with a single 1 = 
amplitude and two separate coherent 1 = 1 amplitudes. 
This is the result we show here. More complete details of 
the fits will be given in the separate paper [47] , but here 
we present the "best fit" results. 

The fits were made to a reduced data set in order to 
exactly match the kaon angular coverage of the three de- 
cay modes, and to remove data points in the vicinity of 
the A(1520) where there was evidence (Fig. 17) of less 
than perfect Monte Carlo matching. There were a total 
of 34 free parameters and 1128 data points. The reduced 
X 2 of the fit was 2.15, the best we achieved with any am- 
plitude combination. Most of the parameters were taken 
up with the overall strength of each amplitude, Cj, in 
each W bin. The centroid, width, and Flatte parame- 
ters of the fitted amplitudes as per Eq. (19) are given in 
Table III. 



TABLE III. Results of the fit using one I 
Breit-Wigner line shapes. 



and two / 



Amplitude 


Centroid 


Width 


Phase 


Flatte 




m R 




A$/ 


Factor 




(MeV/c 2 ) 


(MeV/c 2 ) 


(radians) 


7 


1 = 


1338 ± 10 


85 ± 10 


N/A 


0.91 ±0.20 


1=1 (narrow) 


1413 ± 10 


52 ± 10 


2.0 ±0.2 


0.41 ± 0.20 


1 = 1 (broad) 


1394 ± 20 


149 ± 40 


0.1 ±0.3 


N/A 



The 1 = piece of the reaction was found in the fit to 
be at the threshold. The fit was flexible enough to 
let this centroid move smoothly below threshold if neces- 
sary, but the fit was optimal with the centroid of the 
A(1405), nominally at 1405 MeV/c 2 , pushed down to 
1338 MeV/c 2 . The rising and falling of the line shape 
is controlled by the opening of phase space from thresh- 
old on the low mass side, and by the inflection caused 
by the opening of the NK threshold on the high-mass 
side. The intrinsic width of the 1 = resonance was fit- 
ted to 85 MeV/c 2 . However, we expect this width to be 
poorly determined due to the dominance of the thresh- 
olds above and below the centroid. The Flatte coupling 
parameter is close to unity. A value of 0.91 means there 
is a strong switch-over to the NK decay mode as the 
available energy exceeds this threshold. This switch-over 
is consistent with theoretical expectations [9] . 

Fig. 20 shows only the E + 7r _ data and the correspond- 
ing fit, including the underlying separate isospin curves. 
The black solid curve shows the dominant 1 = line 
shape which is the same for all W bins. It exhibits a 
distinct edge and change in curvature at the NK mass 
fithresh due to operation of the Flatte effect. It is evi- 
dent that the data demand this sort of slope discontinu- 
ity in the £77 distributions. The fit has some problems 
for W = 2.1, 2.2 and 2.3 GeV, where the prominent nar- 
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rowing around 1400 MeV/c 2 is not reproduced. We have 
been unable to find a combination of fit parameters and 
amplitudes that would improve this situation. 

One sees in this and the next figures the "narrow" / = 
1 contribution (dotted lines) plus a second quite "wide" 
contribution (dashed lines). Only the narrow line was 
allowed to have a Flatte break at the NK threshold, but 
not the very wide contribution. The centroid, width, and 
Flatte parameter for this and the other curves are given 
in Table III. 

Analogous to Fig. 20, Fig. 21 shows only the E°7r° 
data and corresponding fit, including the underlying sep- 
arate isospin curves. The 1 = line shape (solid black) 
is three times the E°7r curves (solid blue), as given in 
Eq. (14), apart from the incoherent background. Here 
the effect of using two 1=1 amplitudes can be consid- 
ered. This channel is all 1 = 0, but in accommodating 
the global fit to all channels, the position, strength, and 
width of the single 1 = piece is affected. The fit is less 
good than when fitting the E°7r° final state alone, and 
of about equal qualitative goodness as when using two 
1 = amplitudes and one 1=1 amplitude [47]. The 
E°7r° channel did not help us discriminate which ampli- 
tude combination is superior. 

Figure 22 shows only the E~7r + data with the corre- 
sponding fit (solid green), including the underlying sep- 
arate isospin contributions. In this case the fits are uni- 
formly good across all values of W. The black curves are 
the same in each panel except for their fitted magnitudes, 
which are the same in Figs. 20, 21, and 22 at each W. 

The fit comfortably accommodates Breit-Wigner-like 
1 = 1 structures centered near 1394 and 1413 MeV/c 2 . 
There are no standard quark-model £ states that would 
fit this description. The observation at least tentatively 
suggests evidence for the 1=1, J p = 1/2" , £* state pre- 
dicted in some extensions of the basic quark model [18]. 
However, our fit is a phenomenological parametrization 
of the 1=1 amplitude, and not a direct identification of 
resonant states. 

The broad 1 = 1 structure is hard to interpret because 
it is so wide. It could result from a non-resonant coher- 
ent 3-body amplitude present in the reaction mechanism. 
The fit is substantially better when including this second 
1 = 1 amplitude, in fact, it is crucial for providing the 
separation between the mass distributions in the thresh- 
old region of the three charge states. 

The component curves for the one 1 = and two 1=1 
amplitudes contributing to da /dm are the same in shape, 
but differ in magnitude, on each panel. It is evident that 
the 1 = strength is the largest contribution to the reac- 
tion, but the two 1 = 1 contributions are far from small 
in comparison. The magnitudes of the isospin compo- 
nents as a function of W are shown in Fig. 23. These 
are the real coefficients as per Eq. (19) that enter each 
Breit-Wigner amplitude (in magnitude). Above 2.2 GeV 
the 1=1 strengths combined are as large as half of the 
1 = strength. The relative phase angle of the broad 
1 = 1 amplitude is close to zero with respect to the 



1 = amplitude. This means there is no interference 
between them apart from the Breit-Wigner phase depen- 
dence. However the two 1=1 amplitudes have a large 
phase with respect to each other, as given in Table III, 
and for this we have is no simple explanation. 

We think the work discussed above makes the case that 
the A(1405), as seen experimentally in photoproduction 
on the proton, is not an isolated 1 = resonance centered 
near 1405 MeV/c 2 . The observed line shape (or mass dis- 
tribution) differs in each of the three E7T decay modes, 
which shows that there is substantial 1 = 1 strength in 
the system. We found it necessary to carefully consider 
the opening of the NK decay mode. We have interpreted 
the 1=1 strength in terms of two Breit-Wigner reso- 
nances that interfere with the pure 1 = state A(1405). 
After this was done, we arrived at a satisfactory repre- 
sentation of the experimental results. Even our best fit 
does not reproduce the data fully, and it is difficult to tell 
whether the remaining discrepancies indicate unresolved 
systematic issues with the data or additional physics con- 
tent that we have not identified. 

According to our best fit, a narrow 1=1 amplitude 
is a substantial piece of the overall production strength 
of what has loosely been called the "A(1405)". A wide 
contribution also appears to be needed. The extra 1=1 
strength must have J p = l/2~ in order to interfere as 
it does with the 1 = amplitude, the true A(1405). It 
must be emphasized that this 1 = 1 strength has nothing 
to do with the standard £°(1385) J p = 3/2+ since that 
state was carefully excluded much earlier in the analysis 
process, both by explicit subtraction and by recognition 
that it cannot interfere in the present angle-integrated 
spectra. Although our angular coverage of the hyperon 
decays is not complete, a majority of the range has been 
measured. 

Assuming we are correct in the identification and as- 
signment of quantum numbers of the Breit-Wigner ampli- 
tudes we see, we can discuss them in light of recent theo- 
retical models. First, the low mass of the 1 = amplitude 
is consistent with predictions of a two-pole structure for 
the A(1405), wherein the lower of the two poles is more 
likely to couple to the Y<ir final state. For example, in the 
chiral-unitary model of Ramos, Oset and Bennhold [48], 
the lower-mass pole is at 1390 + «66 MeV/c 2 . However, 
the same analysis predicts a £ with l/2~ at 1579 + z274 
MeV/c 2 , which is not consistent with the structure we 
see. In the model of Oiler and Meissner [5] the 1 = 
lower-mass pole is on two Riemann sheets at 1379 — i28 
and 1433 — ill MeV/c 2 , whereas the 1=1 pole is at 
1444 - i69 and 1419 + z42 MeV/c 2 . Hence this latter 
model is somewhat closer to our results. The meson- 
exchange model of Haidenbauer et al. [12] also predicts 
a two-pole structure for the A(1405) with positions at 
1334.3 + i62.3 MeV/c 2 and 1435.8 + z25.6 MeV/c 2 . The 
lower of these is at the E7r threshold, as found by us in the 
present fit. When we fit with two 1 = line shapes [47], 
however, the higher mass centroid does not match the 
predicted pole position. The resonance pole positions 
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FIG. 20. (Color online) Data and fits for £ + 7r~, with each panel showing a different value of W. Data and fitted shapes are in 
red. The isospin contributions are 1 = (solid black), narrow 1 = 1 (dotted black), and wide 1 = 1 (dashed black). The black 
curves are the same in all panels except for normalization. The vertical dashed lines show the E-zr thresholds on the left, the 
nominal 1.405 GeV location, and the NK threshold. The incoherent background is shown as a thin dashed line (red). 



in the various models do not correspond directly to the 
centroids of Breit-Wigner mass distributions, so these nu- 
merical comparisons are only qualitative. 

To make a further connection to previous theoretical 
work we can make some remarks about previous efforts 
to identify a S*(l/2") state near 1380 MeV/c 2 . As men- 
tioned in the Introduction, positing such a state was a 
consequence of examining several open issues in hadron 
structure using a 5-quark baryon ansatz [18]. In that 
class of models, the dominant configuration of some ex- 
cited baryons consists of two diquarks and an anti-quark 
in a mutual L = or L = 1 state. This in turn can lead 
to low mass, negative parity, isovector states such as the 
one under discussion here. The results we obtained here 
may relate to the observation that the line shape of the 
£(1385) does not conform to its expected P-wave char- 
acter, as discussed in Section V A. If there is indeed an 
admixture of an I = 1 amplitude with J p = 1/2 - at 
nearly the same mass, one can in principle have interfer- 
ences that modify the line shape of the experimentally- 
seen A-7T final state. We have not pursued this question 



further at this time. We emphasize once again that a P- 
wave decay cannot be biasing our results for the iS-wave 
£7r data because we integrate over the hyperon decay 
angles, canceling any interference. 



The CLAS results for the £-7r mass distributions in the 
vicinity of 1405 MeV/c 2 are compelling in the follow- 
ing sense. The mass distribution differences between the 
charge states are large and systematic across our mea- 
sured kinematic space. The need for 1=1 strength 
is inescapable. Furthermore, we have shown that the 
£(1385)(3/2 + ) is not a player in this phenomenology, 
and we have taken care to show that the K* produc- 
tion background also does not play a role. Finally, the 
line shape fits that we have made show that the 1 = 1 
strength is described at least in part by Breit-Wigner 
1 = 1 amplitudes with the masses and widths given in 
Table III. 
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FIG. 21. (Color online) Data and fits for E°7r°, with each panel showing a different value of W. Data and fitted shapes are in 
blue. All the other lines and curves are exactly the same as in Fig. 20. The 1=1 curves are not included here. The incoherent 
background is shown as a thin dashed line (blue). 



X. CONCLUSIONS 

The mass distributions or line shapes of the invariant 
E-7T mass have been measured in the region of the A(1405) 
using CLAS at Jefferson Lab. All three charge combina- 
tions were measured, and the main qualitative conclusion 
is that they are significantly different from each other 
and none is well represented by a simple Breit-Wigner 
line shape. We have shown that the background from 
the £(1385) — > ^ir^ states is small and well-controlled 
by scaling the dominant £(1385) — > An decay. We 
have shown that the interference with K * £ final states is 
unimportant in the sense that the line shape results are 
unaffected. 

Interference of / = and 1=1 isospin channels ap- 
pears to lie at the root of the differing line shapes for the 
three £7r final states. That is to say, there is J = 1 
strength present with the same J p = 1/2 - quantum 
numbers as the A(1405). Amplitude-level fits suggest 
that there may be a £-like state in this mass range, and 
that the centroid of the 1 = A(1405) state lies essen- 
tially at the £7r threshold. This places the A(1405) far 



from the nominal PDG mass value, in a place where £-7r 
threshold effects will have to be understood quantita- 
tively to obtain an accurate picture of this state. From 
the same analysis, even the £°7r° channel, which is purely 
1 = 0, cannot be represented by a relativistic Breit- 
Wigner line shape alone. We find that a channel-coupling 
to the unmeasured NK final state via a Flatte-style uni- 
tarization can lead to a satisfactory shape, and that in- 
deed this channel-coupling dominates the observed mass 
distribution. Thus, we find some signature effect for a 
two-pole picture of the 1 = A(1405), in which the reac- 
tion amplitude couples significantly to both final states. 
However, we see also how the 1 = 1 amplitude adds one 
more layer of complexity to the experimental picture by 
influencing the charged final states. The choice of one 
1 = and two 1=1 amplitudes presented in this paper 
led to the best fit among several choices. Similar results 
were obtained using two 1 = and one 1=1 amplitude, 
which may correspond more closely to current theoretical 
ideas, but these are described elsewhere [47]. 

In addition to the results shown in this paper, the pho- 
toproduction differential cross sections of the A(1405), 
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FIG. 22. (Color online) Data and fits for E~7r + , with each panel showing a different value of W. Data and fitted shapes are in 
green. All the other lines and curves are exactly the same as in Fig. 20. The incoherent background is shown as a thin dashed 
line (green). 
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A(1520), and £(1385), will be presented in a separate 
paper [32]. Also, the same data have been used to di- 
rectly measure the spin and parity of the A(1405), and 
this result will also be presented separately [49]. 

Clearly, both more theoretical modeling of the present 
results and additional experimental data are needed. 
The present work has, we conclude, provided detailed 
line shape results, to which a parametrization with a 
set of Brcit-Wigner amplitudes shows the importance of 
1=1, J p = l/2~ strength centered near 1394 and 1413 
MeV/c 2 , with a dominant 1 = piece very near the Utt 
threshold. 
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